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ABSTRACT 


Context. The X-ray emission from flares on cool (i.e. spectral-type F-M) stars is indicative of very energetic, transient phenomena, 
associated with energy release via magnetic reconnection. 

Aims. We present a uniform, large-scale survey of X-ray flare emission. The XMM-Newton Serendipitous Source Catalogue and its 
associated data products provide an excellent basis for a comprehensive and sensitive survey of stellar flares - both from targeted 
active stars and from those observed serendipitously in the half-degree diameter field-of-view of each observation. 

Methods. The 2XMM Catalogue and the associated time-series (‘light-curve’) data products have been used as the basis for a survey 
of X-ray flares from cool stars in the Hipparcos Tycho-2 catalogue. In addition, we have generated and analysed spectrally-resolved 
(i.e. hardness-ratio). X-ray light-curves. Where available, we have compared XMM OM UV/optical data with the X-ray light-curves. 
Results. Our sample contains ~ 130 flares with well-observed profiles; they originate from ~ 70 stars. The flares range in duration 
from ~ 10^ to ~ 10‘'s, have peak X-ray fluxes from ~ 10“'^ to ~ 10“*' ergcm“^ s“', peak X-ray luminosities from ~ 10^^ to ~ 10**^ 
ergs“', and X-ray energy output from ~ 10**^ to ~ 10^^ erg. Most of the ~ 30 serendipitously-observed stars have little previously 
reported information. The hardness-ratio plots clearly illustrate the spectral (and hence inferred temperature) variations characteristic 
of many flares, and provide an easily accessible overview of the data. We present flare frequency distributions from both target and 
serendipitous observations. The latter provide an unbiased (with respect to stellar activity) study of flare energetics; in addition, 
they allow us to predict numbers of stellar flares that may be detected in future X-ray wide-field surveys. The serendipitous sample 
demonstrates the need for care when calculating flaring rates, especially when normalising the number of flares to a total exposure 
time, where it is important to consider both the stars seen to flare and those from which variability was not detected (i.e. measured as 
non-variable), since in our survey, the latter outnumber the former by more than a factor ten. 

Key words. X-rays: stars - stars: flare - stars: activity stars: coronae - surveys - catalogs 


1. Introduction 

The X-ray emission from flares on late-type (i.e. spectral-type 
F-M, or ‘cool’) stars is indicative of very energetic, transient 
phenomena (lasting typically from minutes to hours), associ¬ 
ated with energy release via magnetic reconnection (see e.g. re¬ 
views by Haisch, Strong & Rodono 1991; Favata 2002; Favata 
& Micela 2003; Giidel 2004; Giidel & Naze 2009). The 2XMM 
serendipitous source catalogue (Watson et al. 2009, hereafter 
‘2XMM’) and its associated database of source-specific data 
products (including spectra and time-series) provide an excel¬ 
lent starting point for a comprehensive and uniform survey of 
stellar flares, with an unprecedented combination of sensitiv¬ 
ity, time coverage and sky coverage. In addition, simultaneous, 
time-resolved UV and optical-band images are available from 


Send offprint requests to: J. R Rye 

* Based on observations obtained with XMM-Newton, an ESA sci¬ 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA. 

** Tables C.l and C.2 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-tA/. Figures C. 1 and C.2 are 
only available in electronic form via http://www.edpsciences.org . 


the XMM Optical Monitor (OM) telescope for a subset of the 
sources. 

In this paper we present a survey of physical parameters (rise 
time, fall time, luminosity, energy etc) of flare events identi¬ 
fied in the 2XMM time-series (otherwise referred to as ’light- 
curves’) of late-type stars in the (Hipparcos) Tycho catalogue (or 
more specifically, Tycho-2, Hpg et al. 2000, hereafter ‘Tycho’). 
We distinguish between those stars observed intentionally as tar¬ 
gets of the observations and those that lie serendipitously in the 
XMM field of view (FOV); the latter provide, in principle, an 
unbiased sample of stellar flares. 

The Tycho catalogue provides a large sample of stars, and 
has been used extensively in conjunction with previous X-ray 
studies, from the ROSAT mission, e.g. for galactic structure in¬ 
vestigations (e.g. Guillout et al. 1998a,b, 1999) and for follow¬ 
up of the detailed properties of active stars (e.g. Frasca et al. 
2006; Klutsch et al. 2008). 

Our catalogue of 2XMM-Tycho flares provides a basis for 
various studies, including: the identification of high-activity 
stars, the search for ‘super-flares’ (c.f. Schaefer et al. 2000), 
statistics of stellar flaring, comparison with large solar flares, 
identification of suitable datasets for more detailed analysis (e.g. 
time-resolved spectra), inputs to modelling and theoretical stud¬ 
ies of stellar magnetic activity and flare mechanisms, and impli- 


1 




J. P. Pye et al.: XMM-Tycho stellar flares 


cations for the stellar-system environment including exoplanets. 
Some of these aspects are discussed further in Sect. 7. 

Section 2 of this paper outlines the XMM observations and 
describes the matching with Tycho stars. Section 3 describes the 
selection of the ‘flare’ lightcurves and the estimation of flare- 
event parameters. Section 4 describes the associated OM light- 
curves available for a sub-set of the observations. Section 5 
presents the results, in terms of flare parameters and distribu¬ 
tions, together with estimates of visibility/detection thresholds 
and correction factors that need to be taken into account in de¬ 
riving statistical properties of the flare sample. Section 6 presents 
X-ray ‘hardness-ratio’ light-curves, useful for identifying spec¬ 
tral changes (and presumed temperature variations) with time, 
during the evolution of the flares. Section 7 discusses the out¬ 
comes of the survey, and compares them to several published 
models and scaling relations. Section 8 summarises the work. 

2. XMM observations and matching with Tycho 
counterparts 

2.1. 2XMM 

Details of the 2XMM catalogue and its construction, and an out¬ 
line of the XMM EPIC X-ray instrumentation and observations, 
are given in 2XMM. We summarise here the key features rele¬ 
vant to the present work. 

The 3491 pointed observations that formed the basis of 
2XMM were distributed relatively arbitrarily over the whole sky, 
with a total coverage of ~ 1%. 38 320 source detections were 
sufficiently bright that 2XMM time-series are available, corre¬ 
sponding to 30498 unique sources, 16% of the total sources in 
the catalogue. Of these detections, 2307 were indicated as vari¬ 
able, from 2001 unique sources, ~ 1% of the total catalogue. 
The median flux (in the ’Total’ energy band 0.2-12 keV) of de¬ 
tections with time-series is ~ 2x10“'^^ erg cm“^ s“'; for the ‘vari¬ 
able’ subset it is ~ 4 x 10“'^ erg cm“^ s“'. As noted in 2XMM, 
the ability to detect variability falls towards lower fluxes. Light- 
curves for each detection have a bin-width that is an integer 
multiple of 10 s (with a minimum of 10 s), and set by the re¬ 
quirement to have an average of > 18 ct/bin for pn and > 5 
for MOS, as computed from the source-detection count rates. 
The median light-curve duration is ~ 25 ks; however, periods 
of high background can reduce the useful time, and the median 
exposure times excluding these high-background intervals are 
~ 12 ks and ~ 16 ks for pn and MOS cameras respectively. 
The variability indicator was based on a simple ;^f^-test against a 
null hypothesis of constancy, with a probability requirement of 
< 10“^, and used only data during ‘Good Time Intervals’ (GTIs, 
see 2XMM Appendix A). This rather stringent threshold was 
chosen so that the expected number of false triggers over the en¬ 
tire set of 2XMM time-series was less than one (2XMM Sect. 8). 

The observational period covered by 2XMM was 2000.09 
to 2007.25, with a corresponding epoch range for the source- 
detection positions in the catalogue (this is relevant for matching 
with high-proper-motion stars, see Sect. 2.3). 

2.2. Tycho-2 

The Tycho-2 catalogue (H0g et al. 2000, hereafter ‘Tycho’) con¬ 
tains positions, proper motions and two-band {Vj, Bj) photom¬ 
etry for the 2.5 million brightest stars in the sky (from observa¬ 
tions with the Hipparcos satellite); it is 90% complete to V ~ 
11.5, but has ~ 50% of its stars in the range 11.5 V 13. The 
stellar surface density on the sky varies from ~ 150 stars deg“^ at 


the galactic equator, to ~ 50 at h = +30° to ~ 25 at the galactic 
poles (H0g et al. 2000). 

The observational period of Hipparcos/Tycho was 1989.85 
to 1993.21, with a mean observational epoch of ~J1991.5; the 
epoch of the Tycho-2 catalogue is J2000.0 (H0g et al. 2000). 

As sources of further information, we have matched the 
Tycho-2 catalogue with the Hipparcos catalogue (ESA 1997; and 
used revised parallaxes from van Leeuwen 2007), the Tycho- 
2 spectral type catalogue (Wright et al. 2003) and photometric 
distances from the N2K project (Ammons et al. 2006). 

2.3. Matching of the 2XMM and Tycho-2 catalogues 

A combined 2XMM-Tycho catalogue was generated by match¬ 
ing (often alternatively referred to as ‘cross-correlating’ or ‘join¬ 
ing’) the sky positions of the objects in the two input catalogues 
using a maximum positional offset of 5 arcsec, a value chosen 
from examination of the 2XMM positional errors (c.f. Watson et 
al. 2009, Sect. 9.5) and from trial matching out to a distance of 
10 arcsec. (The Tycho positional errors are much less than those 
for 2XMM, by more than an order of magnitude.) 

In performing the matching, care was required to perform 
proper-motion corrections to the Tycho positions before com¬ 
parison with 2XMM, in order to account for the (albeit relatively 
small) number of high-proper-motion stars in the 2XMM fields. 
There were three cases where proper motion (as given in Tycho) 
over the time range of 2XMM (~ 7 years) was S 14 arcsec, 
i.e. ^ 2 arcsec/yr, large enough to significantly affect the po¬ 
sition matching even after correction of the Tycho positions to 
the nominal mean epoch of 2XMM; these were: 61 Cyg A (HD 
201091), 61 Cyg B (HD 201092) and HD 95735, each with a 
proper motion of ~ 5 arcsec/yr. Eor these three stars, there was 
an additional complication in that 2XMM was found to have 
multiple ‘unique’ sources corresponding to the proper-motion- 
corrected positions of each of these stars, since clearly the al¬ 
gorithm used in 2XMM to match individual detections to form 
unique sources (2XMM, Sect. 8.1) did not have proper motions 
available. Each star has two corresponding unique-source entries 
in 2XMM1. 

To remove matches of any one 2XMM detection with more 
than one star of a closely-spaced set (all were pairs of stars, with 
a star-star separation of 6 arcsec), the match was considered 
to be (somewhat arbitrarily) with the optically-brighter (lower V 
magnitude) star, while the information regarding the companion 
object was retained for later use if needed. There were 84 detec¬ 
tions, from 49 sources, where this action was taken; 43 of these 
detections, from 23 sources, had associated time-series, with 3 
being flagged as variable. 

Table 1 summarises the results of matching the 2XMM and 
Tycho catalogues. Overall, we define ‘cool stars’ to be those with 
colour B-V > 0.3 (i.e. spectral type EO or later, e.g. Allen 1973); 
these comprise ~ 75% of the totals, comparable with 83% for the 
fraction of stars in the whole Tycho-2 catalogue with B-V > 0.3. 
(Spectral types for individual stars are considered in more detail 
for the 2XMM variable sources, as discussed in Sections 3 and 
7.1.) Hence, there were, on average ~ 1-2 Tycho stars detected 
per 2XMM field, from ~ 10 times that number in total. 


* The 2XMM unique-source reference numbers (SRCID) are: 
177369, 177371 for 61 Cyg A; 177373, 177375 for 61 Cyg B; 91087, 
91088 for HD 95735. The corresponding total numbers of individual 
source detections (DETIDs) for each star are 6, 6 and 2 respectively (ex¬ 
cluding a likely spurious detection [SRCID = 91090] for HD 95735). 
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During the matching process, no account was taken of the 
2XMM source quality information (see 2XMM Sect. 7), though 
this was used later, when examining the individual light-curves 
(see Sect. 3 and Appendix A). 

The main sources of information on the names and astro- 
physical properties of the stars were the SIMBAD^ database and 
the XMM-Newton data products for external catalogue cross¬ 
correlations (2XMM Sect. 6). 

We have estimated the expected number of chance coinci¬ 
dences of Tycho stars and 2XMM sources by creating ‘simu¬ 
lated’ Tycho catalogues having declinations increased by (some¬ 
what arbitrarily) between 72 arcsec and 288 arcsec relative to 
their actual values and matching these with 2XMM at a maxi¬ 
mum positional offset of 5 arcsec. This procedure indicates that 
the false match probability is $ 3% for the entire 2XMM/Tycho 
sample and has a similar value for the subset of (brighter) 
2XMM sources with time-series (and is little affected by restrict¬ 
ing the 2XMM source quality flag SUM_FLAG to < 2). A sim¬ 
ple analytical calculation based on the number of Tycho stars 
in 2XMM fields and the sky-area searched, yielded consistent 
results. 

The total XMM-EPIC viewing time (i.e. observation time 
per star x number of stars) of the 2XMM-Tycho survey is given 
in Table 1, together with values for various subsets. 

2.4. Derived quantities 

We summarise here the calculation of the main additional quan¬ 
tities that were needed for our work, and not contained in the 
original 2XMM or Tycho-2 catalogues. 

Standard, Johnson magnitudes and colours {V, B—V) for the 
stars were computed from the Tycho-2 values {Vt, Bj) using the 
formulae given in ESA (1997; c.f. H0g et al. 2000). 

Distances to the stars were computed, in order of decreasing 
preference, from: (i) Hipparcos trigonometric parallaxes (ESA 
1997, and using revised values from van Leeuwen 2007); (ii) val¬ 
ues gleaned from the literature, for a small number (13) of stars 
of particular relevance to our flare study (individual references 
are noted in Table C.l); (iii) ‘photometric parallaxes’ from the 
N2K project (Ammons 2006); (iv) photometric parallaxes from 
the apparent magnitude V on the assumption of a main-sequence 
(i.e. luminosity class V) object, and utilising the Tycho colour 
(B — V) information to estimate absolute magnitude (My) and 
hence distance modulus. 

2XMM X-ray luminosity (Lx) of each object was computed 
using the derived distance, and Tycho visual luminosity (Lopt) 
directly from My. 


3. Analysis of the X-ray light-curves 

The 2XMM light-curves of all source detections flagged as vari¬ 
able, and matched as above with Tycho stars, were visually 
examined to characterise the variability and to check for po¬ 
tential problems that might give rise to spurious indications of 
variability. Where available, we used the EPIC-pn light-curves, 
due to the generally higher signaknoise relative to MOSl or 
MOS2. The inspection included the background light-curve and 
the plot of time intervals used in the variability test. In addition 
to the light-curves themselves, we also examined the 2XMM im- 


^ http://simbad.u-strasbg.fr/simbad/ 


ages for any evidence of problems such as source confusion^. 
Appendix A summarises the main issues that led to rejection of 
specific light-curves for the purposes of the present work. Of the 
157 ‘variable’ light-curves examined, we removed 29 from fur¬ 
ther consideration (and not all of these were from cool stars). 

Eurther information on all the remaining 2XMM/Tycho stars 
was sought from SIMBAD and from the XMM catalogue- 
crosscorrelation products. We retained for further analysis those 
stars with spectral type E, G, K or M, and those stars lacking 
available spectral types which had B - V ^ 0.3. In addition, 
six earlier-type (B-A) stars^^ lacking detailed information were 
retained for completeness pending further analysis and investi¬ 
gation. At this stage there were 96 variable light-curves from 
76 stars (and 76 unique 2XMM sources). We will refer to this 
subset as the Cool, Variable Sample (CVS). 41 of these stars 
are the intentional target of the XMM observation, leaving 35 as 
serendipitous measurements. 

In order to check the behaviour of the variability test and to 
search for flares marginally below the 2XMM variability thresh¬ 
old, we have visually inspected all 1143 2XMM/Tycho light- 
curves not indicated as variable in 2XMM and for which the 
Tycho B —V > 0.0. There were 19 light-curves found with ap¬ 
parent variability. After detailed investigation of the associated 
background light-curves and GTIs, 7 were not considered further 
due to possible contamination of the source light-curves from 
background variability or confusion with nearby brighter, vari¬ 
able sources. Of the remaining 12 light-curves (from 11 stars), 
two had variability probability ^ 5 x 10'^. The remaining 10 
all showed apparent variability but with much of it outside the 
GTIs. We retained all 12 cases for further analysis; these all had 
B - V > 0.3. We will refer to this subset as the Cool, Low- 
Variable Sample (CLVS). There was some commonality of stars 
between CLVS and CVS. 7 of these stars are the intentional tar¬ 
get of the XMM observation, leaving 4 as serendipitous mea¬ 
surements. 

Table C. 1 summarises the properties of each star and corre¬ 
sponding X-ray source, in the CVSh-CLVS samples. The table 
contains one row per XMM observation. Eigs. 1-5 and 12 - 
17 show examples of the variable X-ray light-curves^, while the 
complete set (CVS-hCLVS, 108 X-ray light-curves) is given in 
Eigs. C.l and C.2. 

3.1. Characterisation of the X-ray variabiiity 

We have characterised the apparent form of the variability in 
each light-curve by visual inspection of the 2XMM time-series 
graphical data product and by inspecting light-curves binned at 
various time integrations to examine individual features such as 
flares, in more detail. The types of variability can be broadly 


^ We used the 2XMM graphical light-curve products and the 
2XMM summary page for each source detection. See e.g. LEDAS: 
http://www.ledas.ac.uk/amie5/amie5.php?action=basic&catname=2xmm 
^ 2XMM SRCID: 26140, 38388, 41960, 130617, 50242, 99128. The 
first four of these displayed clear flare-like events. 

^ Note that the timeseries data products (FITS files and graphical 
light-curve plots) which form part of the 2XMM catalogue dataset are 
not corrected for off-axis angle i.e. for vignetting. Hence, light-curves 
from multiple observations of the same source at different off-axis an¬ 
gles will, in general, need to be vignetting corrected if they are to be 
directly compared (or their time-resolved count rates compared with 
the time-integrated rates [*_RATE values] given in the catalogue). All 
X-ray count rates (and derived fluxes and luminosities) in this paper are 
vignetting corrected to give uniform, on-axis values; this includes the 
X-ray light-curves in the Figures. 
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Table 1. Summary of the 2XMM-Tycho survey characteristics. 


Quantity 

2XMM-Tycho 

All Cool stars 

No. of Tycho stars in 2XMM fields 

~ 26000 

~ 19000 

No. of 2XMM sources matched with Tycho stars 

3042 

2357 

No. of 2XMM detections matched with Tycho stars 

4772 

3499 

No. of 2XMM sources with light-curves 

808 

611 

No. of detections with light-curves 

1393 

933 

No. of X-ray variable sources / stars 

123 / 120 

91/89 

No. of variable light-curves 

157 

118 

No. of variable light-curves after quality checking 

128 

96 a 

No. of X-ray variable stars after quality checking 

85 

76 a 

2XMM summed viewing time (Ms) on Tycho stars: 
for all detected stars (b) 

119 

87 (82) 

for all stars with X-ray light-curves (a) 


29 (24) 

for all stars with X-ray variability (a) 


3.9 (1.8) 

for all stars with flares (a) 


3.0 (1.4) 

for all stars with fully-observed flares (a) 


2.6 (1.2) 

for all stars with flares, S : N > 10 (a) 


2.1 (0.62) 

for all stars with fully-observed flares, 5 : V > 10 (fl) 


1.9 (0.58) 


Notes: Values in parentheses (...) relate to serendipitous observations, i.e. the star was not the target of the XMM observation, (a) CVS sample; 
(b) Includes all relevant 2XMM detections not only those with light-curves. The terms ‘source’ and ‘detection’ are as defined in the 2XMM 

catalogue (Watson et al. 2009, Sect. 8.1). 


placed into the following categories®: flare: corresponding to a 
clear rise above and then fall back towards, a quiescent level; 
trend: a rise or fall in the source count rate over the course of an 
observation; gradual: a rise then fall, or vice versa, but without 
the relatively fast rise usually associated with ‘classic’ flares; in¬ 
determinate: there was no clear form or structure to the variabil¬ 
ity. These definitions are somewhat subjective, and in particular 
it may be noted that: (a) the acceptance criteria for flares will 
need to be accounted for when considering their statistical prop¬ 
erties later in this work; (b) gradual events may in some cases 
be ‘unusual’ flares rather than e.g. due to active-region evolu¬ 
tion or rotational modulation; (c) a downward trend may in some 
cases be the later stages of a flare whose rise and peak have not 
been observed; (d) indeterminate variability may in some cases 
be due to low-level flares which cannot be individually resolved; 
(e) there are some flares where part of the event (one or more of 
the rise, peak or fall phases) has not been observed and hence 
the event cannot be fully characterised. Table 2 summarises the 
results of the variability characterisation; flaring is the dominant 
form of variability, with > 80% of the X-ray variable stars dis¬ 
playing this behaviour. 

We have characterised each X-ray flare using the parame¬ 
ters: count rate^ (above the quiescent level) at maximum (Cmax); 
the corresponding time (fmax); the times (t^, (ff), during the rise 
and fall of the flare, when the count rate is Cmax/e; the quiescent 
count rate (Cq) outside of flaring (usually determined from a pe¬ 
riod close but prior to the start of the flare). Additional, derived 
parameters included the rise- and fall-time, Tr = fmax - tr and 
Tf = ff - fmax- Conversion of count rate to flux used the same fac- 


® We found no examples of periodic variability in our sample. 

’’ These are corrected to on-axis values. If necessary (where the 
original 2XMM timeseries data products were used), an approximate 
correction was made by multiplying the count-rates by the factor 
ca_8_RATE/AVRATE, where ca_8_RATE is the emldetect source count 
rate for the appropriate EPIC camera and AVRATE is the time-averaged 
count rate of the timeseries (recorded in the data-product FITS-file 
metadata). 


tor as the corresponding 2XMM total-band source detection*; X- 
ray luminosity Lx and X-ray emitted energy Ex'^ were then de¬ 
rived as in Sect. 2.4. In order to judge the statistical significance 
of a flare, we computed the signaknoise (S:N) ratio'® above the 
estimated quiescent emission level, over the time interval Ti -h Tf. 
For some of the statistical analyses presented later, we restrict 
the samples to flares which were ‘fully observed’, defined as 
those with measured values for both Ti and Tf, and we refer to 
the ‘duration’ of a flare as Tr + Tf. The parameter values for each 
flare/event are listed in Table C.2. 

4. The XMM optical/ultraviolet light-curves 

Optical and ultraviolet data can yield information about those 
regions of the stellar atmosphere at lower temperatures than are 
seen in X-rays, hence giving a more complete view of the flaring 
process (e.g. Giidel 2004, Sect. 12.14; Mitra-Kraev et al. 2005a). 

For all CVSh-CLVS stars which fell within the FOV of the 
XMM-Newton Optical/UV Monitor telescope (OM) (Mason et 
al. 2001) during the corresponding EPIC observations, we have 
extracted the optical/UV photometry from the OM source lists". 
For many of the stars which were targets of the observations, 
OM. fast-mode data were available, with a basic time resolution 
of 0.5 s, subsequently binned in the pipeline data products to 

10 s, and then further binned at typically ~ 300 - 3000 s 

to match approximately the corresponding EPIC data plots. Eor 
the serendipitous sources, there were more limited data, in the 

® i.e. we used a single conversion factor for each light-curve, namely 
the ratio of the 2XMM total-band count rate and flux for the relevant 
EPIC camera and filter. 

^ Ex was integrated over the time interval t^ -i- Tf, unless otherwise 
stated. 

With S:N defined in the sense S/ Vt, where S is the number of 
counts above the quiescent level, summed over the time interval t, + Tf, 
and T is the total count (i.e. due to flare -l- quiescent emission -l- non¬ 
source background) over the same interval. 

" The OM full FOV is 17 x 17 arcmin^, covering the central portion 
of the EPIC X-ray FOV. During a given observation, coverage of this 
full field may be non-uniform depending on the specification of OM 
observing modes. 
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Table 2. Summary of the X-ray variability characterisation. 


Variability 

type 

Sample 

Number 
of stars 

Number of 
light-curves 

Number of events 

Total Completely 

observed 

Flare 

CVS 

63 (30) 

76 (32) 

133 (39) 

116(34) 


CLVS 

8(4) 

8(4) 

11(4) 

11 (4) 


CVS+CLVS 

70 (34) 

84 (36) 

144 (43) 

127 (38) 

Trend 

CVS 

10(1) 

11 (1) 

11(1) 


Gradual 

CVS 

3(1) 

3(1) 

3(1) 



CLVS 

1(0) 

2(0) 

2(0) 


Indeterminate 

CVS 

6(5) 

6(5) 

6(5) 



CLVS 

2(0) 

2(0) 

2(0) 


Poor bgd sub? (a) 

CVS 

2(0) 

2(0) 

2(0) 


All 

CVS 

76 (36) 

96 (39) 




CLVS 

11 (4) 

12(4) 




CVS+CLVS 

84 (40) 

108 (43) 




Notes: Values in parentheses (...) relate to serendipitous observations, i.e. the star was not the target of the XMM observation, (a) These occur as 
events within multi-event light-curves, and hence do not add to the total number of stars or light-curves. 


form of fluxes integrated over typically ~ 1000 s, represent¬ 
ing individual instrument exposures; we refer to these data as 
imaging-mode or low-time-resolution data. There were often 
more than one waveband filter used sequentially during an ob¬ 
servation, thus complicating the identification of variability in 
the OM light-curves. However, there were ~ 120 light-curves 
(each corresponding to a unique 2XMM DETID) with at least 
two OM photometry data points in the same filter. We restrict 
ourselves here primarily to presentation of low-time-resolution 
light-curves for serendipitous sources, and with the additional 
constraint that the OM data were taken during at least part of the 
X-ray flaring period. This results in five OM light-curves, from 
four stars, as shown in Figs. 1-5; also plotted are the corre¬ 
sponding X-ray data for comparison. 

All available OM data (fast-mode where available, for tar¬ 
get stars only, and low-time-resolution data for the serendipi¬ 
tous sample) are presented in Appendix C.3 alongside the X-ray 
light-curves. 


5. Results 

From this point, we will focus our studies on the set of 116 flares 
from the CVS sample having complete (as far as we can deter¬ 
mine) coverage of the event profiles. Of particular interest are 
the 34 such flares which originated from (30) serendipitously- 
observed stars (Table 2), since we can use these data to estimate, 
or at least usefully constrain, the frequency of flare events at X- 
ray fluxes much lower than previous serendipitous surveys (e.g. 
Pye & McHardy 1983). 

Figs. 6, 7 and 8 summarise respectively the distributions of 
rise time (xr) versus fall time (Tf), quiescent X-ray luminosity 
(Tx,quies) versus peak X-ray luminosity (Lx,peak), and Lx,peak ver¬ 
sus total X-ray energy (Lx) emitted over the time period Tr -H Tf. 
Simulations have been used to estimate errors in the derived pa¬ 
rameters, as noted below. 

The observed flare rise- and fall-times each span two orders 
of magnitude, ranging from ~ 200 s to ~ 20000 s (c.f. Fig. 6). 
The upper and lower boundaries are set by the observational con¬ 
straints: the former by the longest XMM exposure durations, 
the latter by the sensitivity of the instrument and data analy¬ 
sis method in permitting identification of short-duration events. 
Simulations indicate that the uncertainty on and Xf is ~ +20% 



0 10'^ 2x10^ 3x10^ 

Time - MJD 53000.8Q883 (s) 


Fig. 1. (From top) (a) XMM EPIC MOS2 X-ray light-curve for the 
star 2MASS J23163068+7905362 (= IRXS J231628.7+790531), en¬ 
ergy band 0.2-12 keV, 800-s time bins; (b) corresponding XMM OM 
UV light-curve, with filters as noted on the plots; the duration of each 
OM data point is indicated by the horizontal bars. 


for S:N^ 10 and x ^ 2 ks; for shorter x, the uncertainty may 
be larger: ~ +30% for x ~ 1 ks. It is evident that although the 
majority of the flares have a ‘classsic’ profile with x^ < Xf, there 
is apparently a significant number with Xi > Xf (see e.g. Giidel 
2004, Sect. 12.11). However, taking account of the above un¬ 
certainties, low S:N in some cases, and additional uncertainties 
due to some flares being overlapped (‘confused’) in time, leaves 
only one potentially significant case: DETID=46637 (xr ~ 17 ks, 
Xf ss 9 ks), identified with the T Tau-type star V410 Tau, and 
even in that case it might be argued that the apparently long rise 
time could be an artifact of two or more overlapping events. 

The stars’ quiescient X-ray luminosities range from ~ 10^^ 
to ~ 10^^ erg s ', while the flare peak X-ray luminosities range 
from ~ 10^® to ~ 10^^ erg s"' (see Fig. 7), the lower values being 
set by the observational sensitivity. The upper values for Lx.quies 
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10^ 2x10“^ 3x10“^ 4x10”^ 

Time - MJD 53269.12344 (s) 


Fig. 2. As Fig. 1, but for the star HD 14716 observed with PN. 


Detid: 45611 (2MASS J040721 81 -1 21 0033) PNS003 At= 1600 



Fig. 3. As Fig. 1, but for the star 2MASS J04072181-1210033 observed 
with PN and with 1600-s time bins. 


correspond to the maximum typically associated with cool-star 
coronal emission (see e.g. Giidel 2004). The observed flare peak 
X-ray luminosities range from ~ 0.1 to ~ 50 times the quiescent 
levels, i.e. a dynamic range of 500, with a median Tx.peak/Tx.quies 
of ~l-2 depending on the precise choice of sample. Simulations 
indicate that the uncertainty on Lx,peak (or more strictly the as¬ 
sociated maximum count rate) is ±20% for S:N^ 10, with a 
significant bias (measured/true value) ~ 0.7 at S:N~ 10, becom¬ 
ing insignificant (i.e. ~ 1.0) by S:N~ 15. The apparent strong 
correlation shown in Fig. 7 between Lx,peak and Lx,quies, arises 
from a combination of observational bias and intrinsic stellar 
properties (c.f. Audard et al. 2000): weak flares on stars with 
high Lx.quies. though likely to be occurring, are not detectable 
(towards lower right of plot), while very strong flares on low 
Lx.qiiies Stars (towards upper left of plot) have a low occurrence 
rate (though such ‘super flares’ are occasionally seen, e.g. Osten 


Detid: 45610 (2MASS J040721 81 -1 21 0033) M1S001 At= 1600 



Fig. 4. As Fig. 1, but for the star 2MASS J04072181-1210033 observed 
with MOSl and with 1600-s time bins. 


Detid: 156921 (2MASS J1 31 41 1 03-1 620235) M2S011 At= 1600 



Fig. 5. As Fig. 1, but for the star 2MASS J13141103-1620235 observed 
with MOS2 and with 1600-s time bins. 


et al. 2010). Fig. 7 also demonstrates the intrinsically large range 
in Lx,peak at a given Lx,quies. i-S- any given star produces a wide 
range of flare strengths. The distribution of flare luminosities is 
discussed in more detail below. 

The flare total X-ray emitted energy (see Fig. 8) ranges from 
~ 10^' to ~ 10^^ erg. The photon passband of the EPIC light- 
curves (0.2-12 keV) encompasses > 70% of the emitted radia¬ 
tion for coronal emission at the temperatures kT ^ 1 keV gener¬ 
ally measured during stellar flares (and often in quiescence also; 
see e.g. Giidel 2004). (Temperature values of this order are also 
supported by the time-averaged hardness ratios, see e.g. 2XMM 
Sect. 9.7.) As noted later (Sect. 7.4(5)), Fig. 8 indicates an ap¬ 
proximately linear relation between Ex and Lx.peak. implying 
that flare duration is essentially independent of amplitude. 
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Fig. 6. Flare fall-time versus rise-time. Only ‘fully-observed’ flares 
are shown (see text for details). Key to symbols: blue circles: 
Serendipitously-observed stars in the CVS sample; red diagonal 
crosses: Target stars in the CVS sample. 



Quiescent luminosity (erg/s) 


Fig. 7. Flare peak X-ray luminosity versus quiescent X-ray luminosity 
(energy band 0.2 - 12 keV). Symbols as in Fig. 6. 

The maximum distance of the flaring stars in our survey 
is ~ 1 kpc, as illustrated in Fig. 9. This figure also demon¬ 
strates that, as might be expected from the magnitude-limited 
nature of the Tycho catalogue, the maximum distance of the 
stars with 2XMM light-curves (whether detected as variable or 
not) is a strong function of stellar spectral type (as indicated 
by colour B - V). We can also see that for the CVS sample 
the serendipitously-observed stars are, overall, more distant and 
of earlier spectral type (smaller B — V) than the target stars, 
with the median distances being ~ 140 and ~ 35 pc respec¬ 
tively (though with considerable overlap in the distributions); 
the corresponding values for the overall Tycho ‘cool-star’ (i.e. 
B - V > 0.3) light-curves set (irrespective of variability) are 
~ 145 and ~ 45 pc respectively. For all Tycho ‘cool’ stars falling 
in 2XMM fields (irrespective of detection in 2XMM), the me- 



Flare peak luminosity (erg/s) 


Fig. 8. Total X-ray emitted energy versus flare peak X-ray luminosity 
(energy band 0.2 - 12 keV). Symbols as in Fig. 6. 



1 10 100 1000 


Distance (pc) 

Fig. 9. Stellar distance (pc) versus colour (B-V, as an indicator of spec¬ 
tral type). Key to symbols: blue filled circles: Serendipitously-observed 
stars in the CVS sample; red diagonal crosses: Target stars in the CVS 
sample; black dots: all Tycho stars with 2XMM light-curves; grey dots: 
all Tycho stars falling in 2XMM fields (without regard to 2XMM detec¬ 
tion). For the last two classes, some points lie outside the plot and are 
therefore omitted, in order to improve visibility of the main distribution. 

dian distance is ~ 195 pc. For the serendipitously-observed stars 
these figures give some measure of re-assurance that our sample 
of flares is reasonably representative, while for the target stars 
they reflect the obvious fact that the objects were observed as 
known, bright, coronally-active stars. 

As mentioned earlier, the threshold for detection of the flares 
is essentially set by their signaknoise (S:N) relative to the qui¬ 
escent emission level. Hence, the thresholds in terms of peak 
flux, fluence, peak luminosity or emitted energy will vary among 
the datasets depending on the photon-counting noise in the indi¬ 
vidual light-curves. Fig. 10 shows the distribution of flares in 
S:N, with a clear drop in the observed numbers for S:N:^ 8-10. 
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Fig. 10. Distribution of (fully observed) flares in signabnoise (S;N), in 
S:N bins of 2. Grey, solid line: all fully-observed flares (127) in the 
CVS-l-CLVS samples; blue, dashed line; fully-ohserved flares (34) in 
the CVS serendipitous sample. The plotted histograms have been trun¬ 
cated at the right. 

Although it is not possible to give a definitive value for com¬ 
pleteness versus S:N, simulations support the view that flares 
with S:N down to at least 10 should be detected at ~ 100% 
efficiency, and we will use S:N=10.0 as threshold to define a 
‘complete’ sample where necessary, e.g. for estimation of flare 
frequency (Sect. 7.4.2). 

6. X-ray hardness ratios and light-curves 

As additional data products specifically for our flare survey, we 
have generated spectrally-resolved light-curves, utilising X-ray 
‘hardness-ratios’. The hardness-ratio plots clearly illustrate the 
spectral (and hence inferred temperature) variations characteris¬ 
tic of many flares, and provide an easily accessible overview of 
the data, prior to further analysis such as time-resolved spectral 
fitting. 

An X-ray spectral hardness ratio was defined as in 2XMM 
as: HR = (Rj - R ,)/(Rj + Rj), where R, and Rj are the count rates 
in energy bands i and j. Here, we choose as ‘standard’, the bands 
i and j to be 0.2-1.0 and 1.0-12.0 keV, in order to (a) provide a 
good dynamic range of sensitivity to typical temperature ranges 
in stellar coronae (kT ~ 0.1-10 keV), (b) maximise signal-to- 
noise and (c) for consistency with 2XMM (i.e. our bands are the 
sum of two or more 2XMM bands). HR light-curves were gen¬ 
erated specifically for our flare survey. Examples are shown in 
Figs. 12-17. We have generated HR light-curves at various time 
resolutions; the examples presented here are all at a time binning 
of 1600 s, for good signal-to-noise whilst showing the main fea¬ 
tures of the flare events. In Fig. 16(c,d) we present an example 
for a ‘harder’ HR, using the bands 1.0-2.0 and 2.0—4.5 keV, per¬ 
mitting tracking of higher temperatures than our standard band, 
but with reduced signal-to-noise, and hence useful only for rela¬ 
tively bright sources. 

Fig. 11(a) shows predicted HRs, calculated for an optically- 
thin, thermal (coronal) spectrum and negligible line-of-sight 
photoelectric absorption, using the XSPEC software package 
and the MEKAF spectral model (e.g. Dorman & Arnaud 2001; 



Temperature (keV) 



Temperature (keV) 

Fig. 11. (From top) (a) XMM EPIC MOS and PN hardness ratios (HR, 
see text for definition) versus temperature, predicted from an optically- 
thin, thermal spectral (MEKAL) model (using XSPEC), assuming neg¬ 
ligible absorption column density. Solid lines - ‘standard’ HRs; 0.2-1, 
1-12 keV; dashed lines - ‘harder’ HRs: 1-2, 2-4.5 keV. Symbols: MOS 
- square and diamond; PN - triangle and circle. Use of the MOS and 
PN medium-thickness filters was assumed in all cases, but this does not, 
in general, have a major effect on the derived HRs. (b) As (a), but for 
the predicted emission measure (cm^'^pc^^) per unit count rate (ct/s). 

Arnaud, Dorman & Gordon 2010). Fig. 11(b) shows the corre¬ 
sponding count rate to emission measure conversion factors. 

All hardness-ratio timeseries plots are presented in 
Appendix C.4 alongside the X-ray count-rate light-curves. 

7. Discussion 

7. 1. Identifications and stellar data 

Most of the serendipitously-observed stars have little previously 
reported information (as gleaned from searching the SIMBAD 
database and other catalogues available via CDS) other than as¬ 
trometric data and apparent magnitude, colour and a ROSAT 
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Detid: 23485 (HD 14716) Exposure: M1S001 



Time - MJD 53269.12344 (s) 



Time - MJD 53403.73898 (s) 

Detid: 45610 (2MASS J040721 81 -1 210033) Exposure: PNS003 





2x10’ 4x10’ 6x10’ 

Time - MJD 53403.73898 (s) 


Fig. 12 . (From top) (a) XMM EPIC MOSl X-ray light-curve for the star 
HD 14716, energy hand 0.2-12 keV, 1600-s time bins; (b) correspond¬ 
ing hardness-ratio light-curve using bands 0.2-1, 1-12 keV, approximate 
temperatures are indicated, (c), (d) as panels (a, b) but for EPIC PN. 


All-Sky Survey (Voges et al. 1999) X-ray flux. References to 
stellar data are given in Table C.l, where we also cite previously 
published XMM-based work on the target stars. 

In a bid to increase our knowledge of the astrophysical 
properties of the serendipitously-observed stars, we searched 
the WASP (Pollacco et al. 2006) database of visible-band 
light-curves and found five matches with our serendipitously- 
observed stars. Of these, three exhibited no notable variability, 
but two showed a clear periodic signal (West 2009); 2MASS 
J04072181-1210033 (DETID = 45610 and 45611, see Figs. 3, 
4, 13), period 6.1 hr; HD 268974 (DETID=54373, see 
Fig. 14), period ai 4.05 dy. The WASP light-curve shapes 
suggest that both systems may be eclipsing binaries, with 
2MASS J04072181-1210033 being a W UMa-type contact bi¬ 
nary. The latter star has also been independently identified as 
a W UMa system by Farrell et al. 2009. We, and Farrell 2009, 
have also found a; 6.1 hr periodicity in the XMM X-ray data 
of 2MASS J04072181-1210033 (DETID = 45610 and 45611). 
HD 268974 (= ASAS J050527-6743.3 = ASAS J050526- 
6743.3) was found as an eclipsing binary in the ASAS survey 
for periodic variable stars (Pojmanski 1998). 

Other ‘well identified’ objects in the serendipitously- 
observed sample include; 


Fig. 13. As Fig. 12, but for the star 2MASS J04072181-1210033. 


- UY Scl, DETID=1916; W UMa-type contact binary (Chen 
et al. 2006), period 8.6 hr, discussed further below); 

- 78 Tau, DETID=48976; 6 Set-type variable and member of 
the Hyades open star cluster; 

- V807 Tau, DETID=50256; PMS star, Orion-type variable; 

- HD 95735 (=G1 411), DETID=121057; dMe flare star; a 
possible exoplanet (GJ 41 lb) was reported (Gatewood 1996) 
but appears not to have been confirmed; it is not in the ex¬ 
oplanet list at http;//exoplanet.eu/ (Schneider et al. 2011) as 
of 2014 October 7; 

- IM Vir, DETID=148010, 148011; Algol-type eclipsing bi¬ 
nary (Pourbaix et al. 2004), period 1.3 dy; 

- NU UMa (= HD 237944), DETID=118876; eclipsing, spec¬ 
troscopic binary, period 5.5 dy (Otero & Dubovsky 2004; 
Griffin 2009); 

- BN Sgr, DETID=204777; Algol-type eclipsing binary, pe¬ 
riod 2.5 dy (Maikov et al. 2006); 

- EQ Peg, DETID=243985; dMe flare star; note that this ob¬ 
ject appears also in the target-star sample (DETID=243984). 

There are at least 18 stars (i.e. ~two-thirds) of the 
serendipitously-observed sample for which further optical spec¬ 
troscopy and time-resolved photometry would be highly benefi¬ 
cial in determining the detailed nature of the objects. 

We note that, as of 2014 October, there are 10 target stars 
out of a total of 45 in Table C.l for which we can find no pub- 
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Detid: 54373 (HD 268974) Exposure: PNS005 
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Fig. 14. As Fig. 12, but for the star HD 268974. 


Time - MJD 52327.34675 (s) 

Fig. 15. As Fig. 12, but for the star SDSS J033815.04+002926.0. 


lished reports of the XMM results in the refereed literature (from 
a search of the CDS SIMBAD database). 


7.2. Comparison of UV and X-ray light-curves 

Although we cannot draw any broad conclusions from such a 
small and incomplete sample (five UV light-curves from four 
stars. Sect. 4), it can be seen that the UV emission in some 
cases clearly shows flare-like behaviour (e.g. Figs. 1, 2) in oth¬ 
ers the UV flux does not show an obvious correlation with the 
X-ray flare (e.g. Fig. 3) (see e.g. review by Giidel & Naze 2009, 
Sect. 2.4.2, and references therein). Inspection of X-ray and UV 
(fast-mode) light-curves for those target stars where both are 
available, shows a wide variety of relative variability between the 
two wavelength regimes, even after allowing for possible time- 
offsets between X-ray and UV flares of up to several hundred 
seconds (e.g. Mitra-Kraev et al. 2005). 

For the two cases where there was UV coverage both in the 
pre-flare quiescent phase and in proximity to the X-ray flare 
peak, i.e. DETID=242179 (star 2MASS J23163068+7905362 
= IRXS J231628.7+790531, B - V ~ 0.8, Fig. 1) and 
DETID=23485 (star HD 14716, spectral type E5 V, Eig. 2), we 
can, at least tentatively and in a limited manner, apply the anal¬ 
ysis of Mitra-Kraev et al. (2005a) in comparing the X-ray and 
UV luminosity and flare energy. (It may be noted in passing that 
the Mitra-Kraev et al. sample was composed of [five] dMe stars, 


while the two stars from our survey are of earlier spectral types 
[E-G].) The two flares in our survey have large X-ray amplitudes 
(from quiescent to peak luminosity), of a factor ten or greater, 
compared with those discussed in Mitra-Kraev et al. (typically 
amplitudes of a factor 2-3), and are of much longer duration 
(~ 10"^ s versus 10^ s). In contrast, the amplitude of the UV 
flare emission for our two 2XMM-Tycho stars is of order a factor 
two, similar to that in the Mitra-Kraev et al. sample. Mitra-Kraev 
et al. refer to a similar effect, i.e. an apparent deficit in UV emis¬ 
sion or excess in X-rays, for the ‘flat-top’ flares in their sample. 
Eollowing the discussion of Mitra-Kraev et al., comparison of 
the X-ray-to-UV flare-emitted-energy to the X-ray-to-UV quies¬ 
cent (or mean) luminosity for the two 2XMM-Tycho stars may 
suggest for these large-amplitude, long-duration events, the pres¬ 
ence of an additional energy source or mechanism, boosting the 
‘hot’ X-ray-emitting plasma relative to the cooler UV-emitting 
region. As a further investigation into this effect, we have se¬ 
lected a good example of an apparently isolated, large amplitude, 
long duration flare from one of the target stars in our sample; the 
well-known flare star 61 Cyg B (DETID=229736, spectral type 
K7 V), observed in the OM UVW2 Alter (see Sect. C.3). The 
conclusions are the same as those noted above for HD 14716 
and 2MASS J23163068+7905362. 
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0 5000 10"* 1.5x10'' 2x10"' 2.5x10^ 3x10“' 

Time - MJD 53000.80883 (s) 



Time - MJD 52228.05497 (s) 


Fig. 16. (From top) (a), (b) As Fig. 12, but for the star 2MASS 
J23163068+7905362 (= IRXS J231628.7+790531), MOS2 only, 400- 
s time bins, (c), (d) As panels (a, b), but for the ‘harder’ HR using the 
energy bands 1-2, 2-4.5 keV; 1600-s time bins. 


Fig. 17. As Fig. 12, but for the contact-binary, W-UMa-type star UY Scl 
(period 31 ks, Chen et al. (2006)), PN only; two observations, each 
of duration ~ 30 ks, with a gap of ~ 6 ks from the end of the first to the 
start of the second observation. 


7.2.1. The eclipsing RS CVn-type binary SV Cam 

This is one of our target-star sample. Sanz-Forcada et al. 
(2006) have presented the XMM measurements from 2001 
(DETID=75409) and 2003 (DETID=75410), and reported an 
eclipsed X-ray flare during the 2001 observation. They discuss 
only the X-ray data. We now present the corresponding OM UV 
timeseries (see Eig. C.l, DETID=75409) which clearly shows 
the eclipse, and its registration with the X-ray light-curve. An 
eclipse is also seen in the OM UV data for 2003 (see Eig. C.l, 
DETID=75410), though with no obvious X-ray counterpart. 

7.3. Hardness-ratio light-curves 

The HR light-curves demonstrate'^ the ability to track tempera¬ 
ture*^ (and emission measure) changes with time. As indicated 
by the predictions in Eig. 11, and borne-out by the observed flare 


We make the plausible assumption that other physical parameters 
(e.g. absorption column density) that might cause changes in spectral 
shape, are not time variable for these sources. 

Since stellar coronae are, in general, multi-temperature sources, the 
single ‘temperature’ derived from the HR is a ‘mean’ value, weighted 
by the distribution of emission measure with temperature. 


light-curves in Eigs. 12-16, our ‘standard’ HR is sensitive over 
the temperature range kT ~ 0.3-3 keV, whilst the ‘harder’ HR is 
useful over a somewhat higher range of kT ~ 1-10 keV, albeit 
with reduced signal-to-noise. The somewhat different effective 
area as a function of energy for the XMM EPIC MOS and PN in¬ 
struments leads to distinct HR(T) curves in each case, with MOS 
having a rather ‘harder’ response and greater rate of change with 
temperature than PN (c.f. Pig. 11). 

In general, as can be seen in Pigs. 12-16, there is a clear, 
rapid temperature increase at flare onset; subsequent behaviour 
varies from flare-to-flare. In contrast, the flux variability in the 
contact-binary star UY Scl (DETID=1916) in Pig. 17(c,d) shows 
no significant changes in HR. However, the same feature was not 
apparent on the previous rotation of the system (Pig. 17(a,b), 
DETID=1917). Hence it is unclear whether this variability is 
due to flare-like activity or e.g. rapid active-region evolution. 
An analysis of the XMM X-ray spectrum of UY Scl has been 
reported by Stobbart et al. (2006). 
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7.4. Flare parameter distributions and application to flare 
models 

We present here some examples of comparing the 2XMM-Tycho 
results with a selection of previously published flare models (or 
more directly, the associated diagnostic plots), scaling laws and 
observational flare surveys. This is not intended to be exhaustive 
but rather an illustrative comparison. Advantages of the current 
survey over many previously reported ones include uniformity 
of measurements and analysis, and sensitivity. 

The physical parameters of temperature and emission mea¬ 
sure have been derived from the observational quantities: hard¬ 
ness ratio and count rate, using the conversion factors presented 
in Section 6 (Fig. 11). The count-rate to emission-measure 
conversions used were those for the corresponding hardness- 
ratio/temperature . 

We summarise here the main findings. 

1. Temperatures from hardness ratios. These provide a uni¬ 
formly distributed (in time) and relatively high time reso¬ 
lution, useful for identifying and characterising e.g. flare on¬ 
set (e.g. Fig. 12 - Fig. 16) or multiple, overlapping flares 
(e.g. Fig. C.2, DETID=53358,61224,243984). The HRs are 
temperature-sensitive over a limited range, from ~ 0.3 to ~ 3 
keV for the ‘standard’ ones (Fig. 11), and hence the ‘mea¬ 
sured’ temperatures will be restricted, and in particular, will 
tend to underestimate the flare-peak values (since, in general, 
the plasma may be multi-temperature), and hence care must 
be exercised in using them quantitatively. 

We have investigated this effect by comparing our HR- 
derived flare-peak temperatures with those from spectral 
fits reported in the literature for the same XMM observa¬ 
tions. Our comparison set, though rather inhomogeneous 
and incomplete, comprised 19 flares from 12 target stars, 
and clearly confirmed that our HR temperatures consistently 
fall below the spectral-fit values, with the HR temperatures 
~ 0.3-0.5 of the spectral-fit ones. 

2. Emission measures from count rates. The derived emission 
measures (and luminosities and emitted energies) are rel¬ 
atively insensitive to the precise choice of temperature; as 
shown by Fig. 11(b), the conversion factor, for the ‘standard’ 
0.2-12 keV band, varies by at most a factor ~ 2 over the 
range ~1-10 keV, or over the range ~0.2-l keV. 

3. Peak temperature vi peak luminosity, peak emission mea¬ 
sure, and emitted energy. Our results (see Fig. 18(a)) are gen¬ 
erally similar to those of Aschwanden et al. (2008) as sum¬ 
marised in their Fig. 1 (see also Giidel 2004, Sect. 12.12). 
For example, our fltted'^^ power-law slopes of 7.7 and 8.1 
for peak emission measure as a function of temperature, for 
the target and serendipitous samples respectively, are com¬ 
parable with 4.5 + 0.4 of Aschwanden et al., noting that (a) 
the angular difference in inclination for slopes of 4 and 8 
is only 7 deg.; (b) the ‘compression’ of high temperatures, 
previously mentioned, will tend to increase the slope. The 
ranges of temperatures and emission measures covered by 
our sample and Aschwanden et al. are similar, though (again 
probably due to the ‘compression’ effects) the 2XMM-Tycho 
temperatures are, in general, rather lower than those reported 
in Aschwanden et al. We have a larger (by factor ~ 2) and 


We have used the same method as Aschwanden et al. (2008), i.e. 
linear ordinary least-squares bisector to calculate the power-law slope. 
In performing the fits we have discarded some of the extreme, ‘outlier’ 
points. 


more uniform and coherent dataset (i.e. all from 2XMM), 
and clearly confirm that larger flares are hotter (Giidel 2004). 
Our fits ignored ten ‘outliers’ at T > 3 keV; these flares 
have relatively low peak emission measure for the peak tem¬ 
perature, or equivalently, relatively high peak temperature 
for the peak emission measure. They are: 4 flares from the 
classical T Tau-type pre-main-sequence (PMS) stars SU Aur 
(DETID 53358; 3 flares; Robrade et al. 2006, Eranciosini et 
al. 2007b) and CR Cha (120401; Robrade et al. 2006), both 
XMM target objects), and 6 flares from the serendipitous 
sample - HD 31305 (AO type; 53325; Arzner et al. 2007; 
Eranciosini et al. 2007b), TYC 9275-01654-1 (194425), 
TYC 1082-02107-1 (224075), IRXS J231628.7H-790531 
(242179), BN Sgr (E3 V Algol-type; 204777; Maikov et al. 
2006), 2MASS J05350341-0505402 (possibly a PMS star in 
Orion Molecular Cloud 2/3; 64214)'^. An alternative, or pos¬ 
sibly additional, explanation for the high hardness ratios and 
hence high derived temperatures, could be the presence of 
significant line-of-sight absorption, preferentially removing 
low-energy photons from the observed flux. Indeed, Robrade 
et al. 2006 report significant X-ray absorption with column 
density nn ~ 3 x 10^* cm“^ for both SU Aur and CR Cha. 
They also And a hot component with kT ~ 5 keV in SU Aur. 
Eranciosini et al. (2007b) report significant, but lower ab¬ 
sorption, «H ~ 6 X 10^° cm“^, for HD 31305, together with 
a high-temperature flaring component with kT ~ 9 keV. 
Eor the remaining five stars, we have performed spectral 
fits'® to the time-averaged spectra, yielding significant ab¬ 
sorption for BN Sgr (hh ~ 1 x 10^' cm“^) and 2MASS 
J05350341-0505402 (wh ~ 2 x 10^' cm^^), and low absorp¬ 
tion («H ^ 4 X 10^® cm“^) in the other three cases. 

Giidel (2004, Sect. 12.12) has noted that ‘non-flaring’ emis¬ 
sion from a G-star sample of ‘solar-analogues’ follows a 
similar X-ray luminosity / temperature trend to the flare- 
peak properties of his flare dataset, suggesting that flares 
may contribute systematically to the hotter components of 
the plasma. As an extension of this scheme, and somewhat 
speculatively, we have plotted essentially all timeseries data 
in our 2XMM-Tycho variable, cool-star sample in Eig. 18(b), 
where each point represents a single time bin (usually of 
« 400 or » 1600 duration); there are ~ 5500 data points. 
Although there is a wide spread, there appears to be a clear 
trend, similar to that exhibited by the flare-peak values. Eor 
ease of cross-comparison, we also show the flare information 
from Eig. 18(a). 

4. Peak temperature vi duration. Our results are broadly similar 
to those of Aschwanden et al. (2008), who obtain a power- 
law slope of 1.8 + 0.2 with a relatively low degree of correla¬ 
tion (correlation coefficient = 0.39). We And slopes of ~ 2.1 
and ~ 1.0 for our target and serendipitous samples respec¬ 
tively, with correlation coefficients of 0.36 and 0.02. 

5. Peak luminosity vi emitted energy (see Eig. 8). We obtain 
power-law slopes of ~ 1.2, ~ 1.2 and ~ 1.0 for our full sur¬ 
vey, target and serendipitous samples respectively, in good 
agreement with the slope of 1.16 reported by Wolk et al. 
(2005) from their Chandra survey of 41 K-type PMS stars 
in Orion. As Wolk et al. note, this near-linear relationship 
implies that duration is essentially independent of peak lu- 


‘5 The flares from BN Sgr and 2MASS J05350341-0505402 are only 
in the CLVS sample and hence do not appear in Fig. 18(a). 

'* We used the standard XMM data products, fitting with XSPEC and 
a two-temperature MEKAL + uh model. 
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minosity. We have also directly confirmed this in our survey 
dataset by comparing these two parameters. 

6. Time difference between maximum emission measure and 
maximum temperature. At — tEM,„„„ — As discussed by 
Reale (2007), if At > 0, the heat pulse (which originates the 
flare event) is relatively short compared to the characteristic 
cooling time of the emitting loop volume, i.e. it indicates that 
the loop does not reach equilibrium conditions; conversely. 
At - 0 indicates that loop has reached equilibrium. In the 
majority (~ 70%) of cases in our survey we are unable to 
confirm evidence of non-equilibrium (i.e. At is not signifi¬ 
cantly greater than zero), though this may be due to limited 
time resolution and S:N; for ~ 30% there is some evidence, 
but we caution that in most cases the delay corresponds to 
only 1 time bin (ss 400 or a; 1600 s); for ^ 10% the delays 
are several thousand seconds (~ 3-8 ks). (These percentages 
apply to the survey as a whole and individually to the tar¬ 
get and serendipitous samples.) Examples of delays occur 
in the light-curves shown in Figs. 12, 13, 15 and 16. Reale 
(2007) cites 3 example measurements, with At (his Afo-Af) 
~ 20, 0.4, 0.2 ks. 

7. Decay-phase emission measure va temperature. We have 
computed the power-law slope of the emission measure- 
temperature distribution during the flare decay phase, a diag¬ 
nostic of the presence of continued heating during the decay 
phase (e.g. Reale et al. 1997; Reale 2007). However, we have 
not been able to achieve reliable results, probably due to the 
limitations of the single-temperature / hardness-ratio analy¬ 
sis. Hence we defer further consideration, but note that in 
principle, our survey could yield a relatively large sample of 
slope values. 

In summary, the diagnostic plots based on hardness-ratio 
temperatures and associated emission measures etc, provide ex¬ 
cellent, general indicators for the investigation of large numbers 
of flares, but individual, quantitative studies may require more 
detailed analysis e.g. via time-resolved spectral fitting. 

7.4.1. Flare luminosity, energy and emission-measure 
distributions 

It is evident from Fig. 7 that, although flares from the 
serendipitously-observed stars and target stars cover broadly 
similar ranges of both Lx.quies and Lx, peak, the Lx.peak/Lx.quies val¬ 
ues at a given Lx.quies tend to be higher, e.g. ~ 50% of flares from 
the serendipitously-observed stars have Lx,peak/Lx.quies ~ 2, 
while the corresponding fraction for target stars is ^ 10%. This 
may be largely or wholly an observational selection effect aris¬ 
ing from the generally smaller distances and hence higher X-ray 
fluxes of the target stars, enabling the detection of lower levels 
of variability. 

Large, but not infrequent solar flares have soft X-ray peak 
luminosities and emission measures ~ 10^^ erg s“' and ~ 
10^*’ cm“^ respectively, and X-ray emitted energies ~ 10^^ erg 
(see e.g. Giidel 2004, Table 4; Aschwanden et al. 2008, Fig. 3; 
Schrijver et al. 2012), with occurrence rates of ~ 0.1-1/year 
(see e.g. Thomson et al. 2010, Table 1; Schrijver et al. 2012). 
Thus, there is overlap between the lower end of the stellar Lx.peak 
and Lx distributions (Figs. 7, 8) and solar flares. Schrijver et al. 
(2012) conclude that the largest solar flares for which good evi¬ 
dence exists are up to an order of magnitude more energetic than 
those discussed above, placing them well within the 2XMM- 
Tycho distribution. 




Fig. 18. (a) Flare maximum emission measure versus maximum temper¬ 
ature (energy band 0.2 - 12 keV). Symbols as in Fig. 6. (b) As (a) but 
also showing, as black dots, essentially all data points from the 2XMM- 
Tycho light-curves, i.e. each time bin (1 bin = 400 or 1600 s) in each 
light-curve is represented on the plot, irrespective of the presence of a 
flare; the temperature range has been restricted for clarity in showing 
the bulk of the distribution. 


7.4.2. Flare rates and duty cycles 

In a sense, frequency distributions of various flare properties rep¬ 
resent the final distillation of the survey. Measured rates of flar¬ 
ing, or in more general terms, frequency distributions of intrinsic 
physical properties such as peak luminosity and emitted energy, 
allow comparison with models for flare production. The fre¬ 
quency distribution of peak flux is of practical use in predicting 
numbers of flares observable in future surveys and missions. We 
discuss the peak-flux distribution first, together with the topic of 
survey completeness. 


13 






J. P. Pye et al.: XMM-Tycho stellar flares 


Table 3. Summary of the flare frequency statistics and distributions, for the CVS serendipitous-stars sample. 


Quantity 


A 

Raw, Coverage 

measured corrected 

B 

Raw, Coverage 

measured corrected 

Min. selected signabnoise 

10.0 





Min. selected peak flux, /x,peak,min 

(erg cm“’ s“') 

2.5 X 10“*’ 


1.0 X 10“*’ 


Above /x,pcak,min-' 






No. of flares in survey 


15 

16.3 

10 

10.2 

No. of flares predicted, all-sky/year 



4.4 X 10® 


2.8 X 10® 

No. of flares predicted per star/year 

(a) 


4.3 X 10“’ 


2.7 X 10“’ 

Power-law index, Uf 

(b,c) 

0.55 ±0.14 

0.58 ±0.15 

0.94 ± 0.30 

0.94 ± 0.30 

Normalisation, Aref.f 

(b) 

7.0 

7.3 

10.0 

10.2 

with /x,ref = 1 X 10“*’ ergcm“’ s“* 






Flare duty cycle (fraction) 

(a) 


4.3 X 10“’ 


2.0 X 10“’ 


id) 


7.3 X 10“’ 


3.4 X 10“’ 

Min. selected peak luminosity, Lx,peak,min 

(ergs“‘) 

2.0 X 10’® 


1.0 X 10’*’ 


Above ^X,peak,iTiin* 






No. of flares in survey 


14 

16.5 

11 

13.3 

No. of flares predicted per star/year 

(a) 


4.3 X 10“’ 


3.5 X 10“’ 

Power-law index, Ul 

(b,c) 

0.29 ± 0.08 

0.28 ± 0.07 

0.38 ±0.12 

0.37 ±0.11 

Normalisation, Aref,L 

(b) 

00 

bo 

10.5 

11.0 

13.3 

with Lx,ref = 1 X 10’® ergs“* 






Min. selected emitted energy, Fx.min 

(erg) 

1.0 X 10” 


1.0 X 10 ”* 


Above 






No. of flares in survey 


13 

15.4 

9 

11.3 

No. of flares predicted per star/year 

(a) 


4.03 X 10“’ 


2.9 X 10“’ 


(d) 


11.6 


8.5 


(a,e) 


4.9 X 10““* 


4.4 X 10“’ 

Power-law index, 

(b,c) 

0.29 ± 0.08 

0.28 ± 0.08 

0.47 ±0.16 

0.47 ±0.16 


(e) 


0.72 ± 0.20 


0.56 ±0.19 

Normalisation, Aref,E 

(b) 

13.0 

15.4 

26.6 

33.1 

with £x,ref = 1 X 10” erg 







Notes: columns A: results for thresholds (/x,peak,min, Tx, peak,min, £x,min) based on lowest value in the set, with minimum S:N > 10.0; columns B: as 
A, but for substantially higher thresholds, (a) considering all serendipitous-sample stars (504) with 2XMM light-curves, irrespective of detected 
variability; associated survey on-time fmtai ~ 2.4 x 10’s; (b) from ML fit; (c) no correction has been made to these values, e.g. the bias correction 
(N - 1)/N suggested by Crawford et al. (1970), which could reduce the values by up to ~ 10%; (d) considering only serendipitous-sample stars 
(30) observed to flare; associated survey on-time ftotai ~ 1-4 x 10® s; (e) applying scaling factor based on each star’s quiescent X-ray luminosity 
and normalised to a quiescent solar X-ray luminosity, see text Sect. 7.4.2 for details. 


In order to estimate the rate of flaring and the flare ‘duty 
cycle’ we had to take account of possible incompleteness in 
the survey, due to variations in the minimum detectable flare 
strength arising mainly from the source quiescent flux against 
which the flare had to be viewed. Hence, for each flare we com¬ 
puted a ‘survey completeness’ factor, C (range 0.0 - 1.0), be¬ 
ing the fraction of the survey time in which the flare could have 
been detected. Each flare then contributed 1/C to the corrected 
distribution. For the CVS serendipitous sample and our chosen 
S:N threshold of 10.0 (Sect. 5), the maximum correction factor 
1/C was ~ 2 with most flares requiring a correction of 1.3. 
Thus, the corrections applied were relatively small. The calcula¬ 
tion of C is described in Appendix B, and example distributions 
are shown in Fig. 24. 

In order to provide a simple parametrisation of the flare- 
peak frequency distribution, we have considered a (cumulative) 
power-law: 

Ni> /x//x.ref) = A^ref {fx! 

and determined the power-law index Of from the mea¬ 
surements using a maximum-likelihood (ML) method (Jauncey 
1967; Wall & Jenkins 2003). Consistency between measure¬ 
ments and model was checked with the Kolmogorov-Smirnov 
test (e.g. Wall & Jenkins 2003), yielding ZJ^ax ~ 0.25, corre- 

Duty cycle = the sum of the flares’ duration (jr + Tf) as a fraction 
of the sum of the light-curves’ duration (the total ‘on-time’). 


sponding to formal acceptance of the fits at probabilities ~ 20% 
for the serendipitous sample. We have used the power-law form 
purely as a convenient, empirical characterisation rather than 
to demonstrate any fundamental, physically-based shape for the 
measured distributions. 

Table 3 and Figs. 19, 20 summarise the flare-frequency 
statistics, in relation to peak flux (/x.peak) and duration/duty- 
cycle. The values presented are based on only ‘fully-observed’ 
flares; this results in an under-estimate, but this amounts to 
20%. The selected threshold flux /x,peak,min = 2.5 x 10“'^ 
(Table 3, columns A) corresponded to the lowest flux in the sam¬ 
ple; we also show the results for a substantially higher value, 
/x,peak,min = 1.0 X 10“'^ (Table 3, columus B), resulting in an 
increase in the best-fit a but with much larger uncertainties. 
(Increasing /x,peak,min still further would result in a rather small 
number of flares available for fitting, and hence even larger un¬ 
certainties.) Due to the small numbers of stars, and the lack 
of detailed information for most of the serendipitous sample, 
we have not attempted to divide them into different categories; 
hence the distributions and statistics refer to a rather heteroge¬ 
neous collection of stellar types. As the target sample comprises 
a rather arbitrary and ill-defined set of objects (other than all 
being well-established active stars), we have not attempted to 
correct and fit their frequency distributions. 
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-2: model fit to 1 

' ■' 3: corrected, measured 

. 4: model fit to 3, Fmin=2.5E-12 

— • - 5: model fit to 3, Fmin=1.0E-12 


lE-13 lE-12 lE-11 

Flare peak flux, f. (erg cm*^ s*^) 


lE-10 



Flare peak flux, fp {erg cm*^ s*^) 

Fig. 19. (From top) (a) Cumulative frequency distribution of flare peak 
X-ray flux, /x,peak. for the CVS serendipitous-stars sample: black solid 
line - ‘raw’, measured distribution; black dashed line - model distri¬ 
bution from ML power-law fit to raw, measured distribution; red solid 
line with dots - coverage-corrected ‘measured’ distrihution; red dotted 
line - model distrihution from ML power-law fit to corrected distribu¬ 
tion; blue dash-dot line - model distribution from ML power-law fit 
to corrected distribution, with higher minimum flux, /x.min- Each model 
distribution has been derived from the data over the flux range indicated 
by the model (power-law) line, (b) As (a), but for the CVS target-stars 
sample, raw, measured distribution. 


We have examined the sensitivity of the derived values to 
changes in the selection criteria and correction factors for the 
serendipitous sample, as follows. As expected, reducing the 
S:N threshold increased the correction factors to be applied. 
However, S:N thresholds of 8.0, 5.0 resulted in only relatively 
modest increases in corrected flare rates and duty cycles, by 
factors ~ 1.4 and 1.8 respectively. Conversely, increasing the 
S:N threshold to 15.0 yielded a reduction by a factor ~ 0.85. 
The effect of varying S:N threshold on the power-law index was 
also relatively small, with the largest change being ^corrected = 
0.77 + 0.15 for a S:N threshold of 5.0. We have also examined 
the effect of errors in the estimation of the coverage-correction 



Fig. 20. Cumulative frequency distribution of flare duration: blue 
(lower) line - the CVS serendipitous-stars sample; red (upper) line - 
the CVS target-stars sample. The distributions are the ‘raw’, measured 
ones, i.e. not coverage corrected. 


factor C; shifting the C distribution by a factor 2, in the sense 
of worsening the coverage (i.e. C(/x,peak,i) ^ C(2/x,peak,i)), in¬ 
creased the corrected numbers by 20% and increased a by a 
factor ~ 1.3-1.4, but otherwise had little effect. 

We are reporting the activity levels of serendipitously- 
observed stars for a sample defined as being generally ‘active’ 
(i.e. X-ray emission detected, but not necessarily flaring or oth¬ 
erwise variable; Sect. 2.1), and relatively optically bright, and 
incomplete in respect of M-type stars (Sect. 2.2). This is im¬ 
portant to bear in mind if comparing these results to other sam¬ 
ples or surveys. In addition, care must be taken in recognising 
the set of stars to which the various statistics are applicable, as 
indicated by the Notes in Table 3, e.g. the estimated flare fre¬ 
quency when considered only within the set of observed flaring 
stars or only those stars with detected X-ray emission, is obvi¬ 
ously higher than if considered over the whole set of Tycho stars 
falling within the 2XMM survey area. 

We note that the flare duty-cycle values from our survey 
are of the same order, i.e. typically a few percent, as those re¬ 
ported by Walkowicz et al. (2011) from a Kepler white-light 
sample, when considering only those stars observed to flare (and 
hence biasing the result upwards). Our flare-frequency distri¬ 
bution power-law indices ckf are comparable to that of Pye & 
McHardy (1983), who derive a value of 0.8 (+0.4, -0.5) from 
an all-sky survey of fast-transient X-ray sources, at least 60% of 
which were likely cool-star flares (Pye & McHardy 1983; Rao 
&. Vahia 1987). However, we caution against too-detailed com¬ 
parisons given the differing selection criteria and sample types. 

Using the principles outlined above, we have also estimated 
the flare rates in terms of intrinsic properties of the flares: peak 
X-ray luminosity (Lx.peak) and X-ray emitted energy {Ex)- We 
recognise, as noted by Giidel (2004, section 13.5), that there 
are likely to be significant biases and incompleteness in such an 
analysis. In particular, the distributions we present are observed 
distributions, not volume-normalised luminosity and emitted- 
energy functions; the coverage correction C takes explicit ac¬ 
count only of Tycho stars with 2XMM time-series, and we see 
the probable effects of incompleteness in the flattening of the dis¬ 
tributions towards low Lx,peak and Ex due to failure to detect in- 
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Fig. 21. As Fig. 19, but for flare peak X-ray luminosity, Lx,peak- Note that 
this is an obsen’ed luminosity distribution, not a (volume-normalised) 
luminosity function. 


trinsically faint flares towards larger distances (see Figs. 21, 22). 
The results are given in Table 3 and Figs. 21, 22, in a similar way 
to the peak-flux values described earlier. For fix, min = 1.0 x 10^^ 
erg, the power-law index an for the serendipitous sample is 
somewhat below the lower-end values (~ 0.5) in the compila¬ 
tion of Giidel (2004), possibly due to incompleteness towards 
lower Ex- If we raise the serendipitous-sample threshold £x,min 
to 1.0 X 10^"^ erg, the best-fit slope steepens to ~ 0.5 (Table 3, 
columns B). In order to mitigate the possible incompleteness ef¬ 
fects, while recognising the limited number of stars and flares 
available, we also attempted to define an approximately volume- 
limited sample by lowering the S:N threshold to 5.0 and setting a 
maximum distance t/max (c.f. Fig. 24(d)) of 100 pc. This gives 11 
flares with Lx,peak > 10^® erg s“' and 8 flares with Ex > 10^^ erg, 
with an estimated completeness ^ 50% (though the maximum 
applied correction factor 1/C ~ 10), resulting in power-law in¬ 
dices of ffL = 0.52 + 0.16 and qte = 0.75 + 0.26. 

Audard et al. (1999, 2000) have reported EUV flare rates 
from a sample of 10 active cool stars. The estimated rate of flar¬ 


Fig. 22. As Fig. 21, but for flare X-ray emitted energy, E. 


ing for our serendipitous-stars sample, even when restricted to 
the stars observed to flare (i.e. Table 3, note (d)), is much lower, 

by a factor ~ 10-100, than those reported by Audard et 

al. for the seven stars in their sample with measured distribu¬ 
tions extending beyond Ex ~ 10^^ erg. Audard et al. (2000) 
noted a roughly linear relation between flare rate (above a de¬ 
fined Ex threshold) and stellar quiescent coronal luminosity, 
Lx.quies- We do not have sufficient detected flares to yield useful 
frequency distributions for individual stars; however, we have 
formed a ‘scaled’ frequency distribution for the serendipitous- 
stars sample by weighting the contribution of each flare inversely 
according to the quiescent X-ray luminosity of its star, i.e. by 
1/Lx,quies- The resulting oe ~ 0.7 (i.e. Table 3 note (e)), is some¬ 
what steeper than the unweighted value, and still consistent with 
the range reported by Audard et al. (2000) and the literature re¬ 
viewed by Giidel (2004). Schrijver et al. (2012) have compared 
solar flare frequency distributions with those for the five G-K- 
type stars from Audard et al. (2000), scaling the latter (in fre¬ 
quency) by Lo,x,quies/L*,x,quies, using a nominal solar quiescient 
X-ray luminosity Lox.quies = 4.3 x 10^^ ergs^'. They show that 
even with the scaling, the stellar rates are substantially higher. 
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by a factor ~ 100, than the solar distribution. For Fx.min ~ 10^^ 
erg, their scaled rates are ~ 10 events/star/year, compared with 
~ 0.01 solar events/year^®. The corresponding scaled rate for our 
serendipitous sample was ~ 5 x 10^“* events/star/year, somewhat 
below the nominal power-law indicated by Schrijver et al. (their 
figure 3), but arguably within their overall error limits. In sum¬ 
mary, our serendipitous-sample results can span the solar rate 
estimates, depending on the set of stars (and associated total on- 
time) within which the detected flares are considered, i.e. nor¬ 
malising to only those stars (30, ~ 1.4 x 10^ s) observed to flare 
will obviously yield a higher rate than normalising to all stars 
(504, ~ 2.4 X 10^ s) with 2XMM light-curves irrespective of de¬ 
tected variability. In this context, we now consider in more detail 
those stars which did not distinguish themselves by detected flar¬ 
ing within the observations, and compare their general properties 
with those that were observed to flare. 

In Fig. 23 we show frequency distributions of the serendip¬ 
itous variable sample (SV), the serendipitous non-variable sam¬ 
ple (SNV) and the target variable sample (TV), for various prop¬ 
erties of the stars and their associated observations. Visual in¬ 
spection indicates; 

- comparing SV and SNV: no substantial difference for any of 
the plotted quantities, i.e. fxifv, Tx.quies, distance, on-time, 
V, B - V, minimum detectable Ex- 

- comparing SV and TV; TV stars tend to have higher fxifv 
and Lx.quies, smaller distance, and be optically redder (i.e. 
later spectral type). All these differences may plausibly be 
attributed to selection effects in the TV sample. 

The ‘survey coverage’ curves in Fig. 24 show that we would 
expect to detect at least 50% of all flares with Ex ^ 10^^ erg 
and duration lO'* s, i.e. the survey incompleteness within these 
criteria is expected to be no more than a factor two. 

Why was there apparently such a low fraction of the 
serendipitous stars seen to flare? We suggest two obvious, and 
not mutually exclusive, possibilities, (i) The stars observed to 
flare are broadly representative of all the serendipitous-sample 
stars with 2XMM light-curves, and the non-detection of flares 
from > 90% of the sample simply reflects the true (rather low) 
flare rate, (ii) The stars observed to flare have some ‘activity’ 
property which manifests itself in flaring, but not in the general 
properties such as quiescent X-ray luminosity. Such additional 
sources of activity might, for example, arise in magnetic-field 
or tidal interactions between the components of a close binary 
system, or interaction between coronal magnetic field and cir- 
cumstellar material in a PMS star (see e.g. Giidel 2004). We 
emphasise that we are not implying that the stars in the non¬ 
variable sample are intrinsically without flare activity, only that 
they have not produced detectable flares in the 2XMM survey 
observations, and we note that this result is broadly consistent 
with estimates based on Poisson statistics (e.g. Akopian 2013). 
The data used in the current analysis were insufficient in terms 
of total observing time for each star to fully resolve issue (i), or 
to address the related topic of degree of correlation, if any, be¬ 
tween flares closely-spaced in time. However, these scenarios are 
amenable to test via future work: both detailed spectrometry and 
photometry of individual stars, and examination of the growing 
XMM observational database, the latter now having publicly- 
available more than a factor three more observations than for 
2XMM. A preliminary inspection of the additional light-curves 
now available via the 3XMM database shows that several of the 
stars in the SNV sample do exhibit flares. 

Following Schrijver et al. (2012), we have assumed Fboi/Fx ~ 3-5. 


8. Conclusions 

We have presented a survey of stellar flares from a sample of 
Tycho cool (F-M type) stars observed by XMM-Newton. The 
results allow a uniform visualisation and analysis of the XMM- 
Newton data; we have augmented the standard X-ray data prod¬ 
ucts to include hardness-ratio time-series, and where available, 
have presented the associated XMM OM UV/visible data. The 
survey has enabled recognition of new, coronally-active stars, 
observed serendipitously by XMM-Newton. We have demon¬ 
strated the utility of such uniform, relatively large samples in 
the statistical investigation of the physical properties of stellar 
coronae and activity, and potentially as diagnostics of under¬ 
lying mechanisms. In a wider context, stellar flares may, as in 
our own solar system, have significant influence on exoplane¬ 
tary systems (see e.g. Dartnell 2011, Feigelson 2010, Horvath & 
Galante 2012, Melott & Thomas 2011), and may contribute sig¬ 
nificantly to the ‘hard’ (photon energy > 2 keV) X-ray emission 
from the Galactic ridge (Warwick 2014). In addition, for plan¬ 
ning of future wide-fleld X-ray missions and instruments (e.g. 
Osborne et al. 2013), the estimates of flare rates and duty cycles 
provide a useful guide to the expected contribution of F-K-type 
stars to the overall frequency of X-ray transient events. 
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Appendix A: Quality screening of the 2XMM 
light-curves 

We list here the main features found in the data products that 
were likely to compromise the quality of the light-curves with 
respect to validity of the detection of ‘variability’ and its charac¬ 
terisation. These features were identified mainly from visual in¬ 
spection of the graphical products. Identification of any of these 
features in a dataset led to rejection of the light-curve for the 
present work. 

Type of feature and number of 2XMM-Tycho variable detec¬ 
tions rejected: 

- Optical loading (see 2XMM Appendix A): 11 

- Confusion and contamination of source flux by another (usu¬ 
ally brighter) nearby source: 8 

- Problems in the background determination (e.g. poor 
background subtraction, background contamination from a 
nearby source): 7 

- Exposure correction or GTI problem (leading to zeros or 
NULL-values in some of the source light-curve bins): 2 

- Spurious detection (due to a nearby bright source): 1 

Thus in total, 29 out of the 157 2XMM-Tycho variable detec¬ 
tions were rejected. 

Of these 29, 11 had 2XMM SUM_FLAG> 3, and none had 
SUM_FLAG=2 (see 2XMM Appendix D, Sect. D.5 for descrip¬ 
tion of the SUM_FLAG values). 
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Fig. 23. Differential frequency distributions of the serendipitous variable sample (SV, blue lines), the serendipitous non-variable sample (SNV, 
grey lines) and the target variable sample (TV, red lines), for various properties of the stars and their associated observations. The distributions 
are normalised to unit area. (From top left, across) (a) Log [X-ray to visual-band flux ratio, /x//v]. (b) Log [quiescent X-ray luminosity, Lx.quies 
(ergs^*)]. (c) Distance (pc), (d) On-time (s). (e) Visual apparent magnitude, V. (f) Colour, B - V. The distributions in (d) relate to individual 
detections (i.e. each DETID was counted once), while the others relate to individual stars (i.e. each SRCID was counted once), (a) and (b) use 
values averaged over all 2XMM detections of each source (i.e. derived from 2XMM SC_ parameters). 


Appendix B: Calculation of the ‘completeness 
correction factor’, C 

Here we outline the calculation of the correction factor, C, used 
in Section 7.4.2 to account for variations in the minimum de¬ 
tectable flare strength. These arose mainly from the source qui¬ 
escent flux upon which the flare was superimposed. C is the 
fraction of the survey time in which the flare could have been 
detected. We note that our primary aim was to demonstrate that 
our flare rates and related estimates were insensitive to incom¬ 
pleteness at the levels of accuracy needed for our analysis and 
warranted by the limited numbers in our samples, i.e. to esti¬ 
mate C to better than a factor ~ 2 and to utilise it at relatively 
modest correction values, i.e. 1/C^ 2. We wished to generate an 
estimate of C in a simple and rapid manner, and based on quan¬ 
tities directly available from the 2XMM catalogue, rather than, 
for example, engaging in extensive simulations. 

The minimum flare peak flux detectable above a signal-to- 
noise threshold of SN^in is: 

/x.peak.SNmin ~ ‘^^mintr'f^cat(^on/fdur)^ ^ 

where: Cf is the EPIC total-band flux error in the 2XMM 
catalogue (i.e. column ep_8_flux_err). Ton is the observation du¬ 
ration (‘on-time’), Tdur is the flare duration, and S:N is evaluated 
over the time interval Tdur (Sect. 3.1). 


The assumption here is that ixf cat provides a reasonable rep¬ 
resentation of the error on the source time-series. Comparison 
of the EPIC total-band count-rate error with the count-rate error 
derived directly from the individual camera (pn, MOSl, MOS2) 
light-curves indicates that the latter can be a up to a factor ~ 2 
greater than the former (especially for MOS1, MOS2)As dis¬ 
cussed in Sect. 7.4.2, we have examined the effects of such errors 
in C on the resulting flare rates and distributions. We have also 
verified that using (Xf cat, rather than an error estimate based on 
the quiescent count rate from the light-curve, results in 20% 
change in the error value. 

Strictly, the value of Ton is not precisely determined, since 
it varies between cameras (pn, MOSl, MOS2), and a flare (or 
other variability) could be flagged in any of the active cameras. 
Our solution for the present purpose was to use MOS 1 on-time 
if available, else MOS2, else pn. This introduces an acceptably 
small ‘uncertainty’ in Ton, e-g- ^ 30% variation in derived Ton in 
~ 80% of cases. 

C is given by the normalised, cumulative distribution of ob¬ 
servation on-times, i.e.: 

^(^ ./X,peak,SNmin,j) — 2/=l ^on,i/ 2/-1 ^on,i 


Though the individual camera count-rate errors from the catalogue 
are generally in good agreement with the corresponding light-curve val¬ 
ues, with a mean ratio ~ 1.0 ± 0.2. 


18 















































J. P. Pye et al.: XMM-Tycho stellar flares 


where: n is the total number of observations in the sample be¬ 
ing considered (and will, in general, include observations where 
no flares were detected), and the set of [/x,peak,SNmin,i, i — 1,«] is 
ordered by increasing value. 

C can also be expressed in terms of flare peak luminosity or 
emitted energy, via: 

f^X.peak ~ 47rr/ /x.peak 

Ex ~ 47rr/2/x,peakTdui- 

where: d is the source distance. 

Alternatively, C can be expressed in terms of maximum ob¬ 
servable distance: 

^max — VI X,peak/ (A^fx,peak) 

I^max ~ VIx7(4^t^Cpeakt^durJ. 

Example curves are shown in Fig. 24. 
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Fig.24. (From top left, across) (a) Survey completeness or ‘coverage’, C, in terms of the fraction of the survey (or specific samples) in which 
a flare with peak X-ray flux, > /x,pcak,min (ergcm^^ s“') could be detected, for a flare S:N > 10 and a duration of 10“* s; /x,peak,min scales as [S:N 
threshold] / [ Vduration]. Key: blue line: serendipitous variable sample (SV); grey line: serendipitous non-variable sample (SNV); the line for all 
serendipitous stars with time-series is indistinguishable from that for the SNV sample, (b) As (a), but for peak X-ray luminosity, > Lx.peak.min 
(erg s“‘); Lx,peak,min scales as [S:N threshold] / [ Vduration]. (c) As (a), but for X-ray emitted energy, > Lx.min (erg); Lx.min scales as [S:N threshold] 
X [ Vduration]. (d) As (a), but for maximum distance, < (pc) to which a flare with Lx.peak > 10^® or Ex > 10^^ could be detected; Vnax scales 
as duration®'^^ x •y/Lx,peak/[S : N threshold] or duration^®'^^ x VLx/[S : N threshold] . 
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Appendix C: Tables and figures available in electronic form only 
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C.1. The stars in the 2XMM-Tycho flare survey 


Table C.l. The stars in the 2XMM-Tycho flare survey (1 row per XMM observation). 
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C.2. The flares in the 2XMM-Tycho flare survey 


Table C.2. The flares and other time-variability events in the 2XMM-Tycho flare survey (1 row per event). 
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C.3. X-ray and ultraviolet light-curves 


Fig. C.l. All 108 EPIC X-ray light-curves (top panel of each pair), and corresponding OM ultraviolet data (bottom panel of each pair) where 
available. Each pair of plots is labelled at the top with the 2XMM DETID, the star name, the EPIC camera and exposure number, and the X-ray 
time binning Al (s). The conversion factor for count rates measured in the different EPIC cameras is 1 MOS count/s » 3.2 PN count/s. The X-ray 
data are for the Total energy band (0.2-12 keV); the OM waveband filters are indicated towards the right of the plot, and colour-coded. EPIC X-ray 
‘flux’ units are Total-band count/s (in one of PN, MOSl, MOS2 cameras), while OM flux units are ergcm^^ s“' for imaging-mode data and 
count/s for fast-mode data. The conversion factors for OM count rates (count/s) to flux values are: 5.67x 10“'^ (W2), 2.20x 10“'^ (M2), 4.77x10“'*’ 
(Wl), 1.99X 10“'® (U), 1.24X 10“'® (B), 2.51x 10“'® (V) (XMM-SOC-CAL-TN-0019 http://xmm2.esac.esa.int/docs/documents/CAL-TN-0019.pdf 
Table 18.) The plots are ordered by 2XMM DETID; within each page, DETID increases from top to bottom, then right to left, starting at top right. 
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C.4. X-ray count-rate and hardness-ratio light-curves 


Fig. C.2. All 108 EPIC X-ray light-curves (top panel of each pair), and corresponding hardness-ratio light-curves with approximate temperatures 
indicated (bottom panel of each pair). Each pair of plots is labelled at the top with the 2XMM DETID, the star name, the EPIC camera and exposure 
number, and the X-ray time binning At (s). The conversion factor for count rates measured in the different EPIC cameras is I MOS count/s » 3.2 
PN count/s. The X-ray count rates are for the Total energy band (0.2-12 keV), in one of PN, MOSl, MOS2 cameras. The hardness-ratios use 
bands 0.2-1, 1-12 keV. The plots are ordered by 2XMM DETID; within each page, DETID increases from top to bottom, then right to left, starting 
at top right. 
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